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We analyze a possible superconductivity in the hole-doped system of layered hydrogenized 
graphene by taking into account thermal fluctuations of the order parameter. In particular, we 
demonstrate that in the one-layer case the values of the high mean-field (MF) critical tempera- 
ture T^ F ~ 80 — predicted recently by Savini et al [T], do not alter significantly due to the 
fluctuations, and the Berezinskii-Kosterlitz-Thouless (BKT) critical temperature of the vortex su- 
perconductivity is almost the same as the MF temperature at doping 0.01-0.1. We show that in 
the case of multilayer system, when the coupling between the layers stabilizes the superconducting 
phase in the form of fluxon superconductivity, the critical temperature T c can increase dramatically 
to the values ~ 150K, higher than the corresponding values in cuprates under ambient pressure. 
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Introduction - Graphene and its modifications are con- 
sidered to be among the most promising elements of the 
future technological devices, in particular due to their un- 
usual electronic transport properties (see, e.g., a review 
Ref. Possible superconductivity (SC) in graphene re- 
mains one of the most intriguing questions. While there 
are some experimental evidences of SC in graphite and 
graphite-sulfur composites, 3 5 in the two-dimensional 
(2D) graphene system such a possibility of SC was stud- 
ied only theoretically in frameworks of different models 
(see, e.g., Refs. |6- 13J ) , and remains in principle an 
open question. Since different calculations demonstrate 
that inter-electron coupling in graphene, including the 
phonon-mediated coupling, is not too large the critical 
temperature is also not expected to be too high. Op- 
timistic estimations include increase of T c up to 10K 
in the doped graphene due to a van Hove singularity in 
the electron density of states (DOS). [12, Moreover, this 
value may be significantly reduced due to thermal fluctu- 
ations of the order parameter (OP). [13] Recently, it was 
predicted by using first-principle calculations, [T] that the 
critical temperature can be much higher in graphane, a 
fully hydrogenized version of graphene . Indeed, as the 
calculations of the authors of paper pQ show, the electron- 
phonon interaction in this system is much stronger com- 
paring to graphene, due to a Kohn singularity in the 
phonon spectrum. The estimation of the electron-phonon 
coupling constant A ~ 1.3 led the authors to the conclu- 
sion that the critical temperature can be of order 80K- 
90K at rather small values of the hole doping, 0.01 — 0.1. 
These estimations were based on the BCS expression for 
the MF critical temperature. However, in principle in the 
2D case the MF approximation is valid only in the case 
of highly doped system, when the Fermi energy is much 
larger than the typical phonon (Debye) frequency. [TU [15] 



At lower values of doping the MF critical temperature 
is significantly reduced due to thermal fluctuations of 
the SC OP. [15] In this case the real critical tempera- 
ture, correspond to the BKT temperature temperature. 
Below this temperature the OP (its phases) becomes al- 
gebraically ordered, forming the so called vortex SC state 
(see, e.g., a review |16j). Since, in the one-layer graphane 
these fluctuations may reduce T c , in the case of multilayer 
graphane, the inter-layer coupling which usually leads to 
an increase of T c of the coupled vortex (fluxon) SC phase 
(see, e.g., Ref. [T7] and references therein), the high criti- 
cal temperature may be even raised above ~ 90K . In this 
Letter, we systematically analyze the role of the thermal 
fluctuations in the doping dependencies of the critical 
temperature of the vortex and fluxon SC in the case of 
one- and multilayer systems. We show that in the first 
case the fluctuations do not suppress T c significantly at 
doping lager than 0.01, and in the second case the inter- 
layer coupling may lead to a significant increase of T c to 
the values higher than the cuprate critical temperature 
135K, the maximal T c under normal pressure known so 
far. 

Model and main equations. -The secondary- 
quantization effective SC Hamiltonian for the multilayer 
graphane can be written in the following form: 
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where the first and the second terms correspond to 
the in-plane and out-plane components of the electron 
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kinetic energies, the third term corresponds to the 
energy of electron-phonon interaction and the last term 
is the free phonon Hamiltonian. In Eq. (El), ^\„{t, r) and 
ipi a (t, r) are the Heisenberg creation and annihilation 
operators of electron (hole) with spin a = ±1/2 in the 
Zth layer at point r and Matsubara (imaginary) time 
t; m, /i and t are the in-plane electron (hole) effective 
mass, chemical potential and inter-layer hopping energy; 
tim = t{$i,m-i + fii,m+i) is the nearest-neighbor inter- 
layer hopping matrix; g is the electron-phonon coupling. 
We consider the simplest dispesionless (Einstein) phonon 
model, so in this case the free phonon propagator, which 
enters in the last term in Eq. ([IJ, has the follow- 
ing form in the Matsubara frequency representation: 
D(ifl n ) = — oj 2 /[£1 2 + u 2 ), where uj and fl n — 2mrT 
(T is temperature and n is an integer number) are 
the Einstein phonon and Matsubara frequencies. The 
first frequency can be considered as a weighted phonon 
frequency in the case of phonons with dispersion. The 
possible values of lu are discussed below. Since the 
one-layer graphane is a direct band semiconductor (with 
the gap ~ 3.5eV near the T-point), the spectrum of 
weakly-doped free carriers can be approximated by the 
standard k 2 expression (contrary to the linear Dirac 
fermion spectrum in the cases pure and weakly-doped 
graphene). The effective masses of the doped heavy 
and light holes and electrons can be aproximated by 
™>hh = 0.64m e , mih = 0.22m e and m\ e = m e , corre- 
spondingly (m e is bare electron mass). [HI [19] Similar 
to the electron-doped case, in the case of hole doping 
we shall use one-band approximation with the average 
effective mass — 2/[l/m;^ + 1/rrihh] — 0.32m e . We 
choose the value of the coupling g such that the BCS 
coupling A = mg 2 /4ir is equal to 1.3 in the hole-doped 
case. [I] In the electron case, the corresponding constant 
must be much larger, due to larger effective electron 
mass. The value of the phonon frequency ujq can be 
estimated as ~ 0.015eV^ in order to reproduce the 
average BCS value for the MF critical temperature 
T c = 1.14w exp[-l/A] - 90A" obtained in Ref. [1 . The 
value of the remaining free parameter in Hamiltonian 
|l]), the interlayer hopping f,, can be estimated from the 
values of the Slater-Koster hopping parameters between 
the ss-, sp- and pp-orbitals ~ ah 2 /m e d 2 , where d is 
the inter-atomic (inter-layer) distance and a ~ 1 is an 
orbital-dependent coefficient (see, e.g., Ref. [20]). Using 
a recent van der Waals-Density Functional Theory result 
for the distance between the centers of masses of two 
neighboring layers of graphane d — 4.65A ,[3T] one can 
assume t ~ 0.1 — 0.3eV^. In order to derive the equations 
for the SC OP $i(T,ri,r 2 ) = (^(r, n)^ U (r, r 2 )), T c 
and the chemical potential as functions of the carrier 
density and the carrier-phonon coupling, one can use 
the expression for the thermodynamical potential of 
the system as a functional of <J>;(r, ri, r 2 ). We con- 
sider the case of the s-wave (isotropic) pairing and 



weak thermodynamical fluctuations of OP. The last 
approximation is valid in the case when the temperature 
is not too low, which is definitely correct when T 
is close to T c and the doping is not extremely low. 
Then, one can show that the OP can be approximated 
by $ ; (ri,r 2 ) ~ A exp[z0j(ri, r 2 )/2)], where A is the 
modulus of the OP (superconducting gap), and 9i(r) 
is its phase. This phase is proportional to the sum 
of the phases of the carrier operators tpif{r, r\) and 
ipii{i~, r 2 ) and depends on the center-of- mass coordinate 
of the Cooper pair. The thermal phase fluctuations are 
dominant comparing to the fluctuations of A, therefore 
we shall consider only fluctuations of 9 assuming that 
A is constant (see, e.g., Ref. [E]). We also neglect the 
temporal dependence of the OP (which is important 
in the quantum fluctuation regime at T — > 0). Then, 
the thermodynamic potential has the following form in 
the second order (hydrodynamic) approximation in the 
fluctuations of the order parameter: A,V0,T) = 

^°*(/i,A,T) + J( ft A,r)/2E,/<l 2 rW(r)) 2 + 
J||(^,A,T)/2^Jd 2 r[l-cos(0,(r)-0 ; _ 1 (r)) 2 , where 
fl pot is the expressions for the MF (or BCS) part 
of the thermodynamic potential, and J([i,A,T) and 
J\\(fJ>, A,T) are the in-plane and inter-plane stiffnesses 
of the SC phases. The explicit expressions for these 
parameters can be found, for example, in Ref. |17j . In 
the MF approximation, variation of this potential with 
respect to A and [i (neglecting the SC phases) leads to 
the BCS-like equation for the SC gap and the number of 
particles equation: 
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where E(k.,k\\) = ^/£ 2 (k, fc||) + A 2 is the quasi-particle 
energy in the SC state and £(k,feii) = k 2 /2m + 2t — 
2tcos(dfcii) — //is the free hole (electron) spectrum. 
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is the charge density distribution in the momentum 
space. Solution of the system of Eqs. ^ and Q at A = 
gives one the dependence of the MF critical temperature 
T C MF on the model parameters. This solution as well as 
the solution for the real critical temperature T c which is 
defined by the stiffness parameters J and Jn in the case 
of one- and multi-layer systems will be discussed below. 

One-layer case. -The solution of the system of MF 
equations ([2]), Q leads to the following approximate de- 
pendence of the SC temperature on nf. 
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which becomes exact at ep << u>q and ep < ujq, corre- 
spondingly. Since in the 2D case the Fermi energy is con- 
nected with the particle density by the following relation 
ep = h 2 irnf/m, one can estimate that the critical value 
of the carrier density above which the MF temperature 
reaches the doping- independent BCS value 1.14u;oe~ 1 / A 
from ep = ujq ~ 0.015eV. This gives the critical number 
of doped electrons per atom: N atom — 0.0017TO/j/m e (we 
used the value v ce u = 5. 4 A for the volume of the unit 
cell in 7i/ = 2N atom /v ceil)- Since the effective hole mass 
is smaller than the bare electron mass, one can conclude 
that the BCS regime in the one-layer case takes place at 
doping much less than 0.01, and that estimations of the 
MF temperature obtained for doping from 0.01 to 0.1 in 
Ref. pQ are correct. As it was mentioned above, there 
is no long-range SC order in 2D case, which means that 
the OP cannot be constant except the T=0 case. Similar 
to the XY spin model, there possible an algebraic order 
in the system at some temperature T 2D below which the 
phases of the superconducting order parameter become 
algebraically ordered . This temperature is defined from 
T 2D = 7r/2J. Using the expression for the energy J, one 
can show that this temperature is approximately equal 
to ep/8 at low doping and approaches the BCS value 
T C MF with doping increasing at e F > uj , i.e. the fluc- 
tuations are not important at doping larger than 0.01, 
and Tl D ~ T C MF ~ 90K in this case. This suggests that 
the pseudogap phase (PC) i.e. the finite temperature in- 
terval between T 2D and T C MF can be observed only at 
extremely low doping (see Fig. la). It is interesting to 
compare how a much large PG region can merge above 
T 2D = ep/8 = h 2 irnf/8m in some of the cuprates. In 
this case, it was estimated that the effective hole mass can 
be doping-dependent: below N at0 m = 0.1 it is of order 
15m e , and after this value of doping it suddenly drops to 
5m e and smoothly decreases with doping increasing. [2"2"] 
As it follows from the last equation for T 2D , the critical 
temperature for cuprates grows linearly with doping and 
reaches 86K at doping 0.1, i.e. there is a large finitite 
temperature interval between T c and T C MF (which may 
be associated with maximal T 2F/ at optimal doping) for 
doping values below a rather large value 0.1. The main 
reason for this difference in the PG phases is much larger 
value of the hole effective mass in cuprates, which leads 
to a slower growth of T 2D with doping, and therefore 
T 2D "meets" T^ F at much larger values of doping (we 
do not discuss here the reason why superconductivity in 
cuprates starts at finite doping ~ 0.05, which is a topic of 
a separate extended discussion) . The solution of the sys- 
tem of equations ^ for T^ IF together with approx- 
imate solution T 2 — e F /8 is presented in Fig. la. As it 
follows from this Figure, the MF approximation is valid 
at values of doping beginning from less than 0.01, which 
means that the critical temperature T 2D ~ 90K can be 
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FIG. 1: The doping dependence of the critical temperaturein 
the case of one-layer (a) of the multilayered system at different 
values of the inter-layer stifness Jn (b). The model param- 
eters used in the calculations are: A = 1.3, uiq = 0.015eV\ 
ran = 0.32m e .Here and in the next Figure, we have used an 
approximate parabolic hole bandwidth W = 2.7eV, which 
corresponds to the maximal doping 0.2. 



reached already at doping less than ~ 0.01! It is impor- 
tant to notice thatat large doping the solution of the gap 
equation |2]) with the number of particle constraint |3]) 
leads to a higher value of T^ IF ~ 105.7-ftT, comparing to 
the BCS solution T C MF ~ 1.14w e- 1 / A ~ 90K. The rea- 
son for this that the BCS result, which follows from the 
equation ^ at large doping (when fj, ~ £f), is valid only 
at small values of A. 

Another interesting result which follows from the MF 
solution at T = is the possibility of the crossover from 
the BSC superconductivity-Bose-Einstein condensation 
(BEC), or superfluidity regime in the system. Indeed, 
the chemical potential fj, ~ ep — j ej, | / 2 becomes negative 
at doping when the Fermi energy is lower than half of the 
pair binding energy eh = — 2wo exp(— 2/A), which means 
that the system is in the BEC regime in this case [23] 
(see also Ref. [H]). Though, the estimations from e F = 
|e&|/2 with the parameters used above show that such a 
crossover can take place at an extremely low doping, less 
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than 0.001. 

Multi-layer case.- It is possible to show, that in the 
layered system T c can be much larger than the 2D tem- 
perature T 2D , and can reach values 8T 2D at large inter- 
particle coupling or doping it can be equal 8T C . Indeed, 
as it was shown by Horovitz [21] , in a rather general case 
with the inter-layer stiffness J\ \ , the equation for the crit- 
ical temperature has the following form: 



T r . = 8T: 
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where E c is the loss of the SC condensation energy 
in the volume £qcZ , where £o = ftw^/A is the in- 
plane coherence length (see also Ref. [T7], where the 
density-dependent solution for the cuprates was anal- 
ized, and references therein). This energy can be esti- 
mated as SC condensation energy density l/2A(e,p)A 2 
multiplied by the volume £ 2 d.[25] Since the DOS at 
the Fermi level is equal N(ep) — rrih / (2irh 2 ) , one can 
find E c = (e F /2n)(l A 2 /v ce u) ~ e_p/10.87r (as above, 
we choose v ce u = 5AA 2 ). Substituting this expression 
into Eq.Q, and using the fact that T 2D ~ 90 A' and 
J,, - {a/d) 2 J = (a/d) 2 (2/ir)T 2D ~ O.OOOleV, one can 
show that at doping ~ 0.01 — 0.1 the E c terms give a sig- 
nificant contribution in Eq. Q, which may lead to a very 
large increase of the critical temperature, T c ~ 150A, 
comparing to the 2D case (Fig. lb). Since the result for 
the T c is rather sensitive to the values of the parameters 
E c and Ju , we have studied such a dependence by vary- 
ing their values one order of magnitude below and above 
with respect to the estimated values E c q = 0.04eF (at 
doping 0.1) and J|| = O.OOOleF. The results are pre- 
sented in Fig. 2. As it follows from this Figure, even a 
rather modest estimation of the values of the parameters 
results in a significant, by 20% — 50%, increase of the 
critical temperature comparing to the 2D case. 

Conclusions.- In this paper, we have considered super- 
conducting properties of multilayer graphane by taking 
into account fluctuations of the order parameter. We 
have shown that in the single-layer case the BKT criti- 
cal temperature which corresponds to the vortex SC is 
equal to the MF temperature ~ 100A beginning from 
a rather low values of doping less than ~ 0.01. Thus, 
the high critical temperature due to Kohn anomaly ob- 
tained in the MF approximation for this range of doping 
PQ is not altered significantly due to thermal fluctuations. 
In the case of multilayer graphane we have shown that 
the inter-layer coupling, which results in the fluxon SC, 
may lead to a significant increase of T c , comparing to the 
single-layer case. Namely, we estimate that the critical 
temperature may reach values ~ 150A, which is signifi- 
cantly higher than the maximal temperature under am- 
bientl pressure in cuprates. 
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FIG. 2: The condensation energy dependence of the critical 
temperature of the multilayered system (in units of T^ D = 
105. 7K) in the case of different values of the inter-layer stiff- 
ness. The reference condensation energy and stiffness are E c o 
and Jno are given in the text. The other parameters are given 
in Fig.l 
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